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Abstract 

The hydrokinetic approach, that incorporates hydrodynamic expansion of the systems formed 
in A+A colhsions and their dynamical decouphng, is apphed to restore the initial conditions and 
space-time picture of the matter evolution in central Au+Au collisions at the top RHIC energy. 
The analysis is based on the detailed reproduction of the pion and kaon momentum spectra and 
femtoscopic data in whole interval of the transverse momenta studied by both STAR and PHENIX 
collaborations. The fitting procedure utilizes the two parameters: the maximal energy density at 
supposed thermalization time 1 fm/c and the strength of the pre-thermal flows developed to this 
time. The quark-gluon plasma and hadronic gas is supposed to be in complete local equilibrium 
above the chemical freeze-out temperature T^h = 165 MeV with the equation of states (EoS) at high 
temperatures as in the lattice QCD. Below T^h the EoS in the expanding and gradually decoupling 
fluid depends on the composition of the hadron-resonance gas at each space-time point and accounts 
for decays of resonances into the non-equilibrated medium. A good description of the pion and kaon 
transverse momentum spectra and interferometry radii is reached at both used initial energy density 
profiles motivated by the Glauber and Color Glass Condensate (CGC) models, however, at different 
initial energy densities. The discussion as for the approximate pion and kaon mT^-scaling for the 
interferometry radii is based on a comparison of the emission functions for these particles. 

PACS numbers: 25.75.-q, 25.75.Cj, 25.75.Ld 
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I. INTRODUCTION 



At present the majority of the dynamical models of A+A collisions, which describe the soft 
physics phenomena are based on the Landau's idea jl| of space-time evolution of the thermal 
matter formed in the collisions. This approach implies at once the specific space-time scales 
in the problem of nuclear scattering such as a time of expansion, a volume occupied by the 
fireball, hydrodynamic lengths, etc. The only direct tool to measure these femtoscopic scales 
is the intensity interferometry method, often called now as the femtoscopy. The measured 
scales - the interferometry, or HBT radii - are associated just with the homogeneity lengths 
in the rapidly expanding system created in heavy ion collisions So, the comparison 
with experimental data of the space-time scales characterizing such and such dynamical 
model of the system evolution and particle production should be one of the first task of a 
justification of the model and discrimination between different approaches. Nevertheless, 
such a comparison was often ignored since almost all dynamic models, which pretend to be 
complete and therefore describe the evolution of the matter as well as its gradual decay, 
e.g., the hybrid (hydrodynamic plus UrQMD) models failed to reproduce pion out-, side- 
, long- interferometry radii simultaneously with the hadronic spectra at RHIC. Until now 
it was possible to reach only when some artificial parametrization of freeze-out processes, 
e.g. a sudden freeze-out at a fairly large temperature close to the hadronization one ji], is 
utilized. 

In Refs. js-?! the HydroKinetic Model (HKM) for A+A collisions has been developed. It 
combines the advantages of the hydrodynamic approximation, where possible phase transi- 
tions are encoded in the corresponding equation of state (EoS), and microscopic approach, 
accounting for a non-equilibrated process of the spectra formation due to gradual particle 
liberation. The dynamical decoupling is described by the particle escape probabilities in 
inhomogeneous hydrodynamically expanding systems in the way consistent with the kinetic 
equations in the relaxation time approximation for emission function [5j. The method can 
be applied to match correctly hydrodynamics and UrQMD using as the input the locally 
non-equilibrated distribution functions from the HKM. Then one can match these models at 
space-like hypersurfaces related to the late stage of the evolution, escaping thus the incon- 
sistencies connected with an inapplicability of hadron cascade models at very high densities 
and with the causality 
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The HKM method also allows one to take into account a back reaction of particle emission 
on the hydrodynamic evolution that corresponds an account of the viscous effects at the 
hadronic stage of the evolution It is worth noting that found in the HKM ratio of the 
shear viscosity to the total entropy is less than 1/2 in the space-time region of maximal 
hadronic emission jol. An analysis of the QGP evolution within viscous hydrodynamics is 
also a topical problem since the shear viscosity brings an important effect, an increase of 
transverse flows during the evolution [l^. However, until the viscosity of the QGP as the 
function of the temperature becomes clear, this effect is simpler to take into account in the 
phenomenological way, as it is proposed in what follows. 

In this article we apply the HydroKinetic Model (HKM) to an analysis of the 

femtoscopic measurements at RHIC for central Au+Au collisions at the top energy ^/s = 
200 AGeV. Namely, we analyze pion and kaon transverse momentum spectra and the m-r- 
behavior of the pion and kaon interferometry radii to clarify, in particularly, how these 
observables depend on the initial conditions: Glauber and CGC-like. The basic hydrokinetic 
code, proposed in is modified now to include decays of resonances into the expanding 
hadronic chemically non-equilibrated system and, based on the resulting composition of the 
hadron-resonance gas at each space-time point, to calculate the equation of state (EoS) in a 
vicinity of this point. The obtained local EoS allows one to determine the further evolution 
of the considered fluid elements. In the zone of chemical equilibrium, above the chemical 
freeze-out temperature, the EoS is taken in accordance with the lattice QCD results. 

The paper is organized as the following. Section II is devoted to the initial conditions 
(IC) for thermal evolution of the matter in Au+Au collisions at RHIC. In Section HI we 
discuss the EoS of the matter in equilibrated and chemically non-equilibrated zones. The 
kinetics of the system in the non-equilibrium zone related to system's evolution and decou- 
pling is described in Section IV. The underlying hydrodynamic model for both chemically 
equilibrated and non-equilibrated domains is presented in Section V. Section VI is devoted 
to the results obtained and discussions. The conclusions and outlook are done in Section 
VII. 
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II. INITIAL CONDITIONS FOR HYDRO-EVOLUTION OF THERMAL MATTER 



Our results are all related to the central rapidity slice where we use the boost-invariant 
Bjorken-like initial condition in longitudinal direction. We consider the proper time of 
thermalization of quark-gluon matter as the minimal one discussed in the literature, Tq = 1 
fm/c |ll|. 

A. Pre-thermal flows 

If one starts the hydrodynamic evolution at the "conventional time" =1 fm/c without 
transverse flow - since no pressure is established before thermalization - the resulting radial 
flow will not be developed enough to describe simultaneously the absolute values of pion, 
kaon and proton spectra, as well as the anisotropy of elliptic flow in non-cental collisions. 
To describe the observables one needs to start the hydro-evolution at very small initial 
time, r ~ 0.5 fm/c |l2|, where it is difficult to expect the thermalization. This controversial 



situation is overcome due to the results of Ref. [13| where is shown that the initial transverse 
flows in thermal matter as well as their anisotropy, leading to asymmetry of the transverse 
momentum spectra in non-central collisions, could be developed at the pre-thermal, either 
classical field (Glasma) even more efficiency than in 

the case of very early hydrodynamics. So, the hypotheses of early thermalization at times 
less than 1 fm/c is not necessary: the radial and elliptic flows develop no matter whether 
a pressure already established. The general reason for them is an essential flniteness of 
the system in transverse direction. Then the flows of particle number or energy directed 
outward the system cannot be compensated by the inward directed (from periphery to the 
centre) flows. This difference means the non-zero net flows no matter how the collective 
velocity is deflned: according to Ekkart or to Landau-Lifshitz. The further development and 
exploitation of these results were done in Refs. 16l-ll8|. 



The initial transverse rapidity proflle is supposed to be linear in radius rx'- 



Vt = where = J < r|, >, (1) 

here a is the second fltting parameter. Note that the fltting parameter a should include also 
a positive correction for underestimated resulting transverse flow since in this work we did 
not account in direct way for the viscosity effects [l^ neither at QGP stage nor at hadronic 



one. In formalism of HKM |6| the viscosity effects at hadronic stage are incorporated in 
the mechanisms of the back reaction of particle emission on hydrodynamic evolution which 
we ignore in current calculations. Since the corrections to transverse flows which depend 
on unknown viscosity coefficients are unknown, we use fitting parameter a to describe the 
"additional unknown portions" of flows, caused both factors: by a developing of the pre- 
thermal flows and the viscosity effects in quark-gluon plasma. 



B. Glauber-like initial transverse profile 



A simple Glauber model initialization assumes that the initial energy density in the 
transverse plane is proportional to the participant nucleon density 1^, 

p(b, Xt) 



e(b, xt) 



eo- 



(2) 



with po = p(0, 0) and 



p(b, xt) = (T(xt + b/2)5(xT - b/2) + T(xr - b/2)5(xT + b/2)) 



1 



1 — 0"7VAf 



T(xr) 



A 



(3) 



where A is atomic number, equal to 197 for Au+Au collision, and (Tatat = 51 mb(=5.1 fm^) 
is the nucleon-nucleon cross-section at y^J/vw" = 200 AGeV. The impact parameter b = (6, 0) 
is equal to zero, 6=0, in the considered case of central collision. The parameter eo = e(6 = 
0, xq" = 0) is the maximal energy density at the initial moment of thermalization. The 
thickness T(xr) is expressed through the Woods-Saxon distribution profile: 



T(xr) 



ws 



(4) 



where 



i^ws(x) 



exp 



X 



T 



Ra 



+ 1 



(5) 



Here we use that Ra = l.UA^^^ -OMA^^^^ ^ 6.37 fm, 6 = 0.54 fm. Constant a is obtained 
from normalization condition: 

T{-kt)cPxt = A. (6) 
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One can think that transversal Glauber-hke e-profile has been formed to some initial time 
To ~ 0.1 — 0.3 fm/c (see below) when the system is not thermal yet. However, the form of 
the profile is, practically, not modified to supposed thermalization time Tq ~ 1 fm/c because 
the transverse velocities reached to this time are relatively small. At the same time, the 
absolute values of energy density can change significantly because of the strong longitudinal 
expansion. We use the maximal energy density eo at time = 1 fm/c as the second fitting 
parameter. 



C. Initial conditions motivated by Color Glass Condensate model 



Within CGC effective field theory some important physical properties of the field are 
defined by the parameter Ag = g'^fi where g"^ = Anag and /i^ is dimensionless parameter, 
which is the variance of the Gaussian weight over the color charges p of partons. The value of 
Aso is approximately equal to the saturation scale value, Qg, and for the RHIC energies one 
can use A^o ~ ~ 2 GeV^ (20|. According to the results of Refs. 2l|,|22|, (proper) time Tq ~ 
3/ As is an appropriate scale controlling the formation of gluons with a physically well-defined 
energy. At later times the dynamics of the classical YangU-Mills fields produced in nucleus- 
nucleus collisions can be linearized and approximated by that of a system of weakly coupled 
harmonic oscillators. Then one can compute the field amplitudes squared in momentum 
space and find corresponding distribution for the gluon number 



22 



23] for cylindrically 



homogeneous transverse profile. It has the form at pt < l.SA^ and 77 = ^ In 



t-XL 



0, 



dN 



exp 



eff J 



ai 



9 



cff 



-1 



(7) 



where mcff = a2Aso, T^e = a^As] a2 = 0.0358, 03 = 0.465. The constant ai/g"^ will be 
absorbed into factor eo which is our fitting parameter. 

The dependence of the distribution ([7]) on transverse coordinates is constructed as 



follows 



22|: 



a?(xt) = a: 



2 P(b,XT) 

sO 

Po 



(8) 



where the participant density at a particular position in the transverse plane is defined by 

(ED. 
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To define the initial energy density profile we need the partonic phase-space distribution 
fo{x,p) = dN/d^xd^p. Note, that it is associated with the hypersurfaces t = const. To 
express the phase-space density through the values ^2^^^^^^ defined at -^/t^ — a;| = tq, one 
should take into account that the density of partons with momentum p crossing element 
d^a{x) of this hypersurface is 

dN 



d^p 



= p^da^{x)fQ{x,p) 

da(x) 

fo{x,p)TQPT cosh Od'^xxdr], (9) 



where 6 = y — t], y is rapidity of partons (in momentum space). Therefore 

f {x p)- ^ (10) 
' ToUiT cosh 9 d'^xxd'^prd'qdy 

One can formally get the d^N distribution from ([7]) by multiplying it by 5-function: 

dN 



d^PTd^XTdrjdy 



f{T,s)S{y-v)- (11) 



Such a phase-space distribution, corresponding the CGC asymptotic results 



24| . is widely 



used for a description of the initial state in A+A collisions j25|. However, a presence of 
the delta-function in the phase-space density contradicts evidently to the basic principle 
of the quantum mechanics. Indeed, the classical phase-space density has to follow from 



the quantum mechanical one in some limit. The Wigner function /w(a;,p) 



that is 



the quantum mechanical analog of the classical phase-space density f{x,p), satisfies the 



restriction J f^{x,p)d^pd^x < {27rh)~^ (see e.g. 27|, note that the equality takes place for 
a pure state only), here the normalization condition J fY^{x,p)d^pd^x = 1 is supposed. It 
evidently excludes utilization of the delta-function as factor in the structure of the Wigner 
function. Therefore, in order to escape contradiction with quantum mechanics, an another 
prescription, instead of utilization of delta function, should be used for the longitudinal part 
of distribution f{x,p); it can be, for example, the boost-invariant prescriptions used in Ref. 



17|. Following to this recept the smearing of (5-function at hypersurface Tq in (fTTj) as follows 



dJl = ; ( ^eff ^ . (12) 

d'^prd'^XTdrfdy \cosh {rj — y) J 

In this way we fix the phase-space density (ITOl) . This may correspond to quasi-thermal 
averaged partonic distribution which can be reached at moment Tq due to quantum effects 
(uncertainly principle), different kind of turbulences and Schwinger-like mechanism of pair 



production in the pulse of strong color field. It does not mean that the true thermalization 
which should be supported by a permanent mechanism of partonic interactions is reached 
at To ~ 3/As ~ 3 fm/c. 

As a result we use the following form of boost-invariant phase-space distribution for 
gluons at the initial hypersurface Tq: 

_2 ai(romTCOsh6')-^ 

JO - 9 / -/ \ > {^'^) 

exp ^y m^g(xr) + Pt cosh 9 / T^si^r) 



here 6 = rj — y, :s.t = {X,Y) = [xt cos (p , Xt sin ip) and we consider gluons as massless 
particles, = p^- Such a distribution depends on the effective mass meg (x^) = a2As(x7-) 
and the temperature Teff(xx') = g.^ A,, (xr ) (numerical values for 02 and a-^ are the same as in 



Eq- dZD), which, in accordance with Ref. [22j, are determined by the local scale As(xq") ([H]). 

The components of the energy-momentum tensor in the pseudo-Cartesian coordinates 
reads 

T^'^(x) = J p^p''f{x,p)pTdpTdyd<j), (14) 

where the Lorentz-invariant integration measure d^p/po in the Cartesian variables is already 
re-written in Bjorken variables as pTdpTdyd(j). 

We numerically calculate the components of the energy-momentum tensor with the dis- 
tribution function, following from Eq. (fT3!) . at 77 = 0. 

Note that, at r = tq, the energy- momentum tensor takes the form 

rn^T, XL = o) = -^KH^tV^ (15) 

where t^'^ are the constant coefficients fixed by the constants 02 and 03. Therefore, the 
energy profile in transverse plane at Tq in central collisions can be presented in the form (see 

m) 

<Xt) = 6o ^'^'^3°f^\ (16) 

Po 

where the number of participants is defined by ([3]). Under the same reason as for the 
Glauber-like IC we use the form of this profile to build the IC for hydrokinetic evolution at 
the thermalization time Tj = 1 fm/c. The maximal energy density eo at (proper) time Tj is 
the fitting parameter as in the case of the Glauber IC. 
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III. THE THERMAL MATTER IN A+A COLLISION AND EQUATION OF 
STATE 



Here we describe the matter properties and its thermodynamic characteristics, e.g. equa- 
tion of state, that are necessary components of the hydrokinetic modeL We suppose that 
soon after thermahzation the matter created in A+A colhsion at RHIC energies is in the 
quark gluon plasma (QGP) state. Also at time Tj, there is a peripheral region with relatively 
small initial energy densities: e(r) < 0.5 GeV/fm^. This part of the matter ("corona") does 



not transform into QGP and have no chance to be involved in thermahzation process 



28|. 



By itself the corona gives no essential contribution to the hadron spectra |28|. One should 
consider it separately from the thermal bulk of the matter and should not include in hydro- 
dynamic evolution. Therefore we cut the initial Glauber or CGC-like profiles at e(r) < 0.5 
GeV/fm^ when consider IC for hydrodynamic evolution of the system. 

During the system evolution the QGP is cooling and finally transforms into hadron phase, 
most probably, according to the crossover scenario. Such a transformation may occur in 
the interval of the temperatures 170-190 MeV. At the temperature T = T^h ~ 165 MeV 
;he chemical freeze-out happens, as demonstrates an analysis of the particle number ratios 
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30|. The conception of the chemical freeze-out means that at the temperatures T > Tch 
the bulk of the expanding matter is in the local thermal and chemical equilibrium while 
at T < Tch the chemical composition becomes in some sense frozen: one can neglect the 
majority of inelastic reactions except for decays of resonances and recombination processes. 
The hadronic matter in the later thermodynamic region is not in the chemical equilibrium, 
moreover, the hadronic medium gradually emits particles being in this zone and, so, loose, 
in addition, also the local thermal equilibrium. Therefore, one should consider in different 
ways the matter evolution in the two 4D space-time zones separated by the 3D hypersurface 
corresponding to the isotherm T = T^h ~ 165 MeV. Let us describe the thermodynamic 
properties of matter in both these regions. 

A. The EoS in the equilibrated space-time domain. 

At high temperatures corresponding to the QGP phase and crossover transition to hadron 
phase we use a realistic EoS jsil adjusted to the lattice QCD results for zero barionic chemical 
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potential so that it is matched with an ideal chemically equilibrated multicomponent hadron 
resonance gas at Tc = 175 MeV. To take into account a conservation of the net baryon 
number, electric charge and strangeness in the QGP phase, one has first to make corrections 



to thermodynamic quantities for nonzero chemical potentials. As it is proposed in [32|, a 
modification of the EoS can be evaluated by using of the Taylor series expansion in terms 
of the light and strange quark chemical potentials, or analogously in baryon and strange 
hadronic chemical potentials: 

p(T,/XB,/ig) _ p(T,0,0) Ixb (I^bY IXs ri^sY Xbs I^b fJ^s 
2^4 2^4 2T^VT/ 2T^VT/ T T 

The expansion coefficients xb and xs are the baryon number and strangeness susceptibilities 
which are related to thermal fluctuations of baryon number and strangeness in a thermal 
medium at zero chemical potentials. 

To obtain the EoS in the equilibrium zone we use the numerical results for xb and xs as 



a function of the temperature given in 



32| . The values for the ratios fiq/T in ( ITTj) during the 



system evolution can be determined approximately. If at some hypersurface corresponding 
to an isotherm, like as at the chemical freeze-out hypersurface, the chemical potentials are 
uniform, then the following ratios remain constants 

^ = consign where q = B, S, E 

during the chemically equilibrated isoentropic evolution of the Boltzmann massless gas. In 
our approximation we use there constrains and find the corresponding constants from the 
chemical potentials obtained together with Tch from an analysis of the particle number 
ratios. In concrete calculations we use the chemical freeze-out temperature T^h = 165 MeV, 
corresponding chemical potentials fiB =29 MeV, fis =7 MeV, fiE =-1 MeV and also the 
strangeness suppression factor 75 = 0.935 which are dictated by 200A GeV RHIC particle 
number ratios analysis done in the statistical model 
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3Q|. 



B. The EoS in the chemically non-equilibrated domain. 



At the chemical freeze-out temperature Tch the "lattice" EoS taken from |3l| and cor- 
rected for non-zero chemical potentials is matched with good accuracy with ideal Boltzmann 
hadronic resonance gas which includes = 359 hadron states made of u, d, s quarks with 
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masses up to 2.6 GeV. Essentially, we use the same particle set in the FASTMC event gen- 



erator 33|. Technically, in the numerical code, we input the corresponding N functions - the 
densities of each hadron i and the equations for already at the very beginning of the 
system evolution; however, these densities are meaningless in the QGP phase and their eval- 
uation does not influence on the system evolution in the equilibrated zone. These functions 
are brought into play at T < Tch- If this thermodynamic region would correspond to the 
complete conservation of the particle numbers then, in addition to the energy-momentum 
conservation, one would account for the conservation equations for particle number flows in 
the form: 

d^{n,u^) = 0, t = l...N (18) 

In our problem, however, during the system evolution in the non-equilibrated zone T < 
Tch the resonance decays have to be taken into account. The decay law in a homogeneous 
medium with T <^ mj (mj is the resonances mass) implies a summing up of a decrease of 
unstable ith particle number due to decays and an increase because of decays of heavier jth 
resonance into ith particle: 

^ = -r.iv, + 5^6,,r,iv, (19) 

j 

where is the total width of resonance i, bij = BijMij denote the average number of ith 
particles coming from arbitrary decay of jth resonance, Bij = Tij/Tj^tot is branching ratio, 
Mij is a number of ith particles produced in j ^ i decay channel. The set on N equations 
( IT9l) . solved together, takes into account all possible cascade decays z —t- j —t- A; . This 

also conserves net charges, e.g. baryon, electric charge and strangeness, since the charges 
are conserved in resonance decay process. If one relates the Eq. (fT9|) to the fluid element 
of some volume AV moving with four-velocity u'^, then a covariant relativistic extension of 
the decay law for a hydrodynamic medium leads to the equation f fTSj) : 

a^(ni(x)M^(x)) = -r,n,(x) + J2 ^yrj^j(a;) (20) 

j 

when one neglects a thermal motion of the resonance j, that can be justified because post 
(chemical) freeze-out temperatures are much less than the mass of the lightest known res- 
onance. Also, Eq. (l20l) for the hydrodynamic evolution is written under supposition of an 
instant thermalization of the decay products, that is consistent with the ideal fluid approx- 
imation (mean free path is zero). In the kinetic part of the HKM we consider the next 
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approximation when the non-equihbrium character of the distribution functions and the 
kinetics of resonance decays are taken into account. We also can approximately account 
for a recombination in the processes of resonance decays into expanding medium just by 
utilizing the effective decay width Ti^eff = I'^i in Eq. (1201) . We use 7 = 0.75 |3J| for the 



resonances containing u and d quarks supposing thus that about 30% of such resonances are 
recombining during the evolution. 

The equations fl20|) together with the hydrodynamic equations and the equation of state 
should give one the energy density and composition of the gas in each space-time points. To 
find the EoS p = p{e, {ui}) for the mixture of hadron gases we start with the expressions for 
energy density and particle density for ith component of multicomponent Boltzmann gas : 

^ -mlT{2,TK2{m,/T) + mK^{mi/T)) exp{fii/T) 



= ^mlTK2{m,lT) exp(/i,/r). (21) 



Then, the equation for the temperature is: 

6 = 3nT + ^n,m / '' J 22 

where n = Having solved this equation numerically for given e and {n^}, we get the 

temperature and then find the pressure using simple relation for multicomponent Boltzmann 
gas: 

p = nT (23) 

The equations (122|) . f l23|) define p = p{e, {ui}). 

Thus, we follow the evolution of all N densities of hadron species in hydro calculation, 
and compute EoS dynamically for each chemical composition of N sorts of hadrons in every 
hydrodynamic cell in the system during the evolution. Using this method, we do not limit 
ourselves in chemically frozen or equilibrated evolution, keeping nevertheless thermodynam- 
ically consistent scheme. 

As it was mentioned before, we use the Boltzmann approximation in the EoS calculation 
to decrease computational time. However, for emission function and spectra calculation 
we use quantum Bose-Einstein/Fermi-Dirac distribution functions with chemical potentials 
calculated to give the same particle densities as in the Boltzmann case. We checked that the 
measure of relative divergence in the energy density if one uses the quantum distribution 
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functions instead of the Boltzmann one, is not bigger than 3% in the thermodynamic region 
which is actually contributed to formation of hadronic spectra. 



IV. KINETICS IN THE NON-EQUILIBRIUM HADRONIC ZONE 

To describe the non-equilibrium evolution and decay of hadronic system we start from 
the Boltzmann equations for the mixture of hadrons, most of which have finite lifetimes and 
decay widths compatible with particle masses. The set of such equations for i-components of 
the hadron resonance gas which account for the only binary interactions (elastic scattering) 
and resonance decays are: 

f = ^r** - Lr'\x,p) + Gt'^'ix^p) - Lt^'ix.p) = G,{x,p) - L,{x,p). (24) 

Here we ignore the processes of resonance recombination which simpler to account phe- 
nomenologically (see the previous Section). The term gain (G) describes an income of the 
particles into phase-space point {x,p) due to scatters and resonance decays. The term loss 
(L) is related to a decrease of particles in the vicinity of the phase space point {x,p) due to 
re-scattering and decays of resonances. The loss term is proportional to the particle number 
density in the point x and so L|™**(x,p) = fiRi, Lf^'^"'^ {x ^ p) = fiDi where R is scattering 
rate, and D is decay rate. If one considers the equations for stable or quasi-stable particles, 
thenLf''^^(x,p) = 0(A = 0). 

The method allowing to find the emission function of the hadrons based on the Boltzmann 



equations in the (generalized) relaxation time approximation was proposed in Refs. 
Following to this method we put: Ji{x,p) ^ Ri^i,eq{x,p), Gi ^ Ri,i.eq.{x,p)fi,i.eq{x,p) + 
Gf^"^""^ {x , p) . The quantity R{x,p) = T~^]{x,p) is the inverse relaxation time, or collision 
rates in global reference frame. Then, 

^d^fi{x,p) = {fi^'^-{x,p) - fi{x,p))Ri{x,p) + Gf''"^(x,p) - Lf'="^(x,p) (25) 

The explicit form of Gf^'^'^^ {x,p) term will be derived later. In the first approximation 
to hydro-kinetic evolution the parameters of the local equilibrium distribution function 
fi,i.eq.{XiP)i G-g- the temperature T{x), chemical potentials fii{x) are determined by the 
hydrodynamic evolution. The details of hydrodynamic approach used in the model are 
described in the next section. 
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A. Emission functions in hyperbolic coordinates and spectra formation 



All our results are related to the very central rapidity interval, ?/ ~ 0, and we will use 
the boost-invariant approach to describe strong longitudinal matter expansion observed at 
RHIC. For such an approach the hyperbolic coordinates in (t, xl) directions are more suitable 
than the Cartesian ones. Then the kinetic equations take a form 

1 / d mr sinh & d d \ 

— mrcosh6'- — +Pt^^ /^(r, t*, r^, Pt) = 

niT cosh y \ or t orj ott J 



fl-^'-ir, 9, TT, Pt) - Mr, 9, tt, Pt)J i?.(r, 9, tt, Pt) + Gr"^(r, 9, tt, Pt) (26) 

where r = -^/t^ — x| is a proper time, rrix = \/ rn^ + p\ is a transverse mass, 9 = rj — y, rj 
is a space-time rapidity, defined above Eq. ([7]), and y is a particle rapidity. 

The formal solutions of f l26|) correspond to the non-equilibrium distribution functions in 
expanding and decaying multi-hadronic system: 

MT,9,rr,PT) = fnro,9<^^'>\T),rp^\T),pT)exp I- J /?,(., ^(^)(r), 4^)(r), Pt)^/^ 



decay / 



TO 



/[^'^■(A,^W(r),4")(r),pr)i?,(A,^W(r),r(.')(r),pT) + Gf'^'^^(A,^W 

(27) 

exp f Ri{s,9^'\T),4\T),pT)ds 



X 



here i?,(A, r^, Pt) = Si?.(A, ^, r^, p^), Gf ^'^^(A, ^, r^, p^) = SGf'^'^^(A, ^, r^, p^). 
Here we use the notation 

sinh 6'(^oHr) = ^ sinh 6* 



, , "° , (28) 

rf»)(r) = FT - ^(rcosh^ - ^r^ + sinh' 9) 



rriT 

The invariant value is p^Ri{x,p) = p*^R*{x,p), where the asterisk * denotes a value in 

the local rest frame of the fluid element in point x, so 

~ cosh y coshj/jAv P^^m p.^ 

= = = ^^9"^^^^^"^^ ^29) 

To connect the formal solution (1271) with observables, e.g. particle spectrum, we use the 

equality 

P^=n ^-r= / da,p>'fix,p) (30) 



d^p 2TxpTdpTdy 

(Tout 

14 



where aout is a "distant" hypersurface of large r = const, where all the interactions among 
hadrons are ceased. 

In what follows we use the variable substitution in the first term of (l27|l describing the 
"initial emission" : 



sinh^ = ^sinh^' 



= + £|(r cosh 6 - ^/r^ + sinh^ 9) 



(31) 



and the substitution : 



sinh 9 = - sinh 9' 



= + ^(rcosh^ - + sinh^ 



(32) 



in the second term of (127|) related to the "4-volume emission". After transformation to new 
variables {r, 6'', F} we arrive with the result: 

/ oo N 

da^p>^f{x,p) = J rfaoV/i"^-(ro,^',r^,p)exp I - J i?,(s, ^(^)(ro), rf ^(ro), Pt)^^ 



CTout 
T 



dX j da,{X)p^ [fl-^'^iX, 9\ r^, Pt)^.(A, 9\ r^, p^) + Gf ^'^^(A, 9\ r^, p^) 

TO o-(A) 

■exp i-^j Us,9^^\X),4\X),VT)ds\ =p''^'^ 



(33) 



d^p 

where 0"(A) is r = A = const hypersurface, so dafj_{X)p^ = Xm^ cosh.9'd9'd'^r'j,. The expo- 
nential values in these expressions are the escape probabilities 



/ oo \ 

P(r,rT,e,PT) =exp i- J R,{s,9^'\T),4\T),pT)ds 



(34) 



for particles with momentum p at space-time point (r, r-p, t] = 9 + y) (in hyperbolic coordi- 
nates) to become free without any collision js, 6|. 

In the expression above we can separate the 4-volume emission function 

^,(A,^,rr,PT) = [fnX,9,rT,p)R,iX,9,rT,p) + G^''^'iX,9,rT,PT)] P(A,rr,^,Pr) (35) 
and the initial emission : 



Si,o{0, Tt, Pt) = fi'^'iro, 9, r^, Pt)'P{to, r^, 9, Pt) 



(36) 
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These expression demonstrate obviously that the particle emission is formed by the par- 
ticles which undergo their last interaction or are already free initially. These expressions 
for the hadron emission function are the basic functions for calculations of the single- and 
™u,«- par«c.e spect.a g To evaluate these ,ua„tit,es for observed (quasi) stable parfdes 
one needs to find the term gain resonance decays and the collision rates i?,. 



B. Resonance decays in multi-component gas 

We suppose that in the first (hydrodynamic) approximation the products of resonance 
decays which interact with medium are thermalized and they become free later, after the 
last collision with one of other particles. However, at the late stages of matter evolution the 
system becomes fairly dilute, so that some of these produced particles get a possibility to 
escape without any collisions: V > 0. To describe this we use the following form for lJ^'^'^°-y 



and Q'^'^'^°-y terms (for 2-particle resonance decay) 35|: 

Pk J Pi Fi^ki 
where resonance i decays into particles or resonances k and /. 

= E E / ^ / ^r,^../,(x,j,,)-^5(^)(p, -Pk- Pd (38) 

j f. J Pj J Pk ^j~^ik 

where the resonance j decays into particles i and k with partial width Tj^ik for this decay 
channel, and 

F,^.k = [ ^S^'\P, -Pk- P.) = —,{{m^ - ml - - ^m\mlfl- (39) 

J Pk J Pi 

To escape the complicated problem with satisfaction of thermodynamic identities in 
hadron resonance gas we utilize in what follows the mass shell approximation for resonances, 
supposing that rrii = {rrii). Also, as it was already discussed, we take into account that the 
resonance mass in hadron resonance gas is much larger than the temperature, rrii » Tc. 
Then the most probable velocity of resonance in the rest system of a fluid element is small, 

pa a/ — , and one can use the approximation 

^ rriiu''. (40) 
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So the resonance distribution function takes the form 

^0 



fj{x,pi) ^ ^nj{x)5^{pj - mju{x)), (41) 



rrij 



It allows us to perform integrations in (138|) over pj, pk analytically and get : 

Gt^'{x,p,) = Y.I1 r,^^k^^^^S{m,u'{x) - pI - p°) (42) 

where p^ = ^/rn^ + (mju(x) — Pi)^. 

Just this form of gain term is used when spectra are evaluated according to Eq. (155]) . 
Note that in practical calculations we substitute (5-function by its Gaussian representation: 

and take a finite parameter value R = 50MeV. 



C. Collision rates 

The collision rate R{x,p) = ^ J^^^^ is one of the basic value for calculation of the intensity 
of the interactions in the expanding system and its decoupling. The latter is described 
through the escape probability V{x,p) (134|) - the integral value of R along the possible 
trajectory of a particle with momentum p running freely through the whole expending 
system. The rate of collisions in the rest frame of some fluid element that accounts for 
scatters of given particle with any other ith hadronic species in the thermal Boltzmann 
system depends only on particle energy E* = p^u^ and the thermodjTiamic parameters of 
this fluid element |36|: 



Here gi = (2ji + 1), Ep = a/ p*^ + m^, Ek^i = A/kfH-mf, Si = [p* + kiY is the squared cm. 
energy of the pair, and o"j(s) is the total cross section of selected particle with particle i in the 
corresponding binary collision. One can change the integration variable to squared center 
of mass energy s, energy of scattering partner E^ and momentum angle 0, and perform E^. 



17 



and ^-integration analytically, which gives the expression for remaining integral: 

no 



is — — m^)E* 



p 



X sinh ( — ^\/ ~ ~ "^lY ~ 4m^m2 ] exp 



1 1 p 



(44) 



We calculate ai{s) in a way similar to UrQMD code 37| : 



Breit-Wigner formula is applied for meson-meson and meson-baryon scattering: 
(^tota/(v^) = (jg, mg, jM, muW Jr, Mr) ^'^ vtc^ ^ i ^ 

^ R^MB^total 

X ■ 



where T total = ^ ^r^mb, with A/i-dependent parametrization of partial decay 

(channels) 

widths: 



M \PcMs{Mr)J 1,09 f PCMsiM)_Y 

^ ^ \vcms{Mr) ) 

chosen to depend on absolute value of particle momentum in two-particle rest frame: 



1 



In the case of meson-meson scattering a constant elastic cross section of 5 mb is added 
in order to fully reproduce the measured cross section. 

• PDG table data for p — p, p — n, p — p, etc. scattering 

• other baryon-baryon scattering: additive quark model: 

= « (I)"""" (l - 4^) (l - 4^) l^bl . 

rrii = 1(0) corresponds to meson (baryon), Si - number of strange quarks in hadron i. 

Note that all relevant resonance states (see above), 359 different species - are taken into 
account for the calculation of crj(s). 
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V. HYDRODYNAMICS 



We describe the system evolution in the equihbrium zone at T > T^h by the perfect 
hydrodynamics. The small shear viscosity effects, which lead to an increase of the transverse 
flows |lO| we account phenomenologically including this effect in the parameter a of initial 
velocity as described in Section IIA. The matter evolution in this zone is described by the 
relativistic hydrodynamical equations related to the conservation of energy-momentum: 

d^T^"" = (45) 

and equations associated with the net baryon number, strangeness and isospin conservations 

d.iqiu") = (46) 

Here qi is the density of conserved quantum number. 

At T < Tch the equations for the system evolution in the first approximation: fi = f\'^'^ 
can be derived from the basic equation (125|) . Namely, integrating the left and right hand 
sides of Eq.f l25p over d^p one arrives to the equation (120]) for particle number flow in the 
non-equilibrium zone, and also to hydrodynamic equation fH5]) by integrating Eq. f|25p over 
Pid^Pi and summing over index i. 

Note that in the first approximation the matter evolution is described by the equations 
of ideal hydrodynamics while the distribution function (127|) in decaying system is non- 
equilibrium. The iteration procedure, including the next order approximations, that, in 
fact, leads to viscous hydrodynamic evolution, is described in Ref. jo]. In this article we 
limit ourself by the first approximation. Then the energy-momentum tensor T^^'^ has a simple 
structure which is employed in this model : 

T^^ = (e + p)uV - p ■ g''\ (47) 

where e is energy density and p is pressure defined from the EoS. In the chemically equi- 
librated zone the pressure is defined from the lattice QCD calculations as discussed in 
Section IIIA. In the non-equilibrium zone the EoS generally depends on all 360 variables, 
p = p(e, {rij}), and it is evaluated altogether with solution of the evolutionary equations. 
The reason is that it is impossible to store the EoS table, therefore we compute pressure each 
time we need it (e.g. when restoring thermodynamic variables from conservative variables 
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or when computing fluxes through each cell boundary) solving analytic equations (122|) . ( 123|) 
numerically. 

Let us rewrite equations in hyperbolic coordinates. These coordinates are suitable for 
dynamical description at RHIC energies, since, for example, zero longitudinal flows corre- 
spond to boost-invariant expansion (so nonzero longitudinal flow correspond to deviation 



from boost invariance), and evolution parameter, r = \/t'^ — is not affected by strong 
longitudinal flow, which saves computational time. It is convenient to write the equations 
in conservative form, then the conservative variables are : 



Q 



I Qr \ 

Qx 

Qy 

Qr) 
\{Qn,}J 



( 7^(6 -hp) -p\ 
V {7^.} J 



(48) 



Here the expression in curly brackets denote variables associated with the particle densities 
for each sort of hadrons. The Qi are conservative variables in the sense that integral (discrete 
sum over all cells) of Qi gives the total energy, momentum and particle numbers, which 
are conserved up to the fluxes on the grid boundaries. The velocities in this expression 
are deflned in LCMS (longitudinally co- moving system), and related to velocities in the 
laboratory frame as : 

cosh Hf 



lab 



,lab 



cosh(yj — 77) 

cosh yf 
cosh(?/j — 77) 



(49) 



where yf = |ln[(l -|- f^"*)/(l — v'-^^)] is the longitudinal rapidity of fluid element, rj 
I ln[(t + z)/(t — z)] is the geometrical rapidity. 
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The full hydrodynamical equations are : 




( Qr \ 

Qx 

Qy 
\{QnAj 



dxP 




v 



+ 



( (g, + p)(i + u2)/r\ 

Qx/r 
Qy/r 

\ {Qnjr} j 



(50) 



fluxes 



and V = [d^, dy, ^9^). 

For hydrodynamic calculations related to midrapidity region on central A+A collisions 
we impose longitudinal symmetry and cylindrical symmetry in transverse direction. This 
actually means that tangential (in transverse direction) and longitudinal velocities in LCMS 
vanish, so Q^j, = Q?, = 0, as well as the fluxes in (p and t] directions. Then, one has to solve 
the following set of equations : 

^ {Qr +P){1 + V^^)/t - {Qr + p)Vr/r ^ 
Qr/r - QrVr/r 
{QnJr - Qn,Vr/r} 




( {Qr+p)Vr^ 



\ 



QrVr + P 
{Qn^Vr} 



I 



\ 



(51) 



fluxes 



Practically Vrjr is ambiguous at r = 0, so we put Vrjr = a there and use a value interpolated 
from the neighboring points. 

We base our calculations on the finite-volume approach: we discretize the system on 
a fixed grid in the calculational frame and interpret as average value over some space 
interval i, which is called a cell {i is a multi-index in multidimensional case). We also split 
continuous time evolution into a sequence of finite timesteps n. 

The Q^ are then updated after each timestep according to the fluxes on the cell interface 
during the timestep Atn- In 3-dimensional case one has the following update formula : 



Qijk~ 



At 



Axi 



i+l/2,jk+Fi_i/2jk 



At 
A^2 



(-^i,i+l/2,fc+-^i,j-l/2,; 



At 



Ax- 



ij,k+l/2+Fij^k_l/2 



(52) 



Where F is the average flux over the cell boundary, indexes +1/2 and —1/2 correspond to 
right and left cell boundary in each direction. 

This is a base of the Godunov method 39| that implies that the distributions of variables 



on the grid are piecewise, this forms the Riemann problem at each cell interface. Then the 
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flux through each cell interface depends only on the solution of single Riemann problem, 
supposing that the waves from the neighbouring discontinuities do not intersect. The latter 



38| 



is satisfied with Courant-Friedrichs-Lewy (CFL) condition [40 1. 

To solve the Riemann problems at each cell interface we use relativistic HLLE solver 
which approximates the wave profile in the Riemann problem by single intermediate state 
between two shock waves propagating away from the initial discontinuity. Together with the 
shock wave velocity estimate, in this approximation one can obtain analytical dependence 
of flux on initial conditions for Riemann problem, which makes algorithm to be explicit. 

We proceed then to construct higher-order numerical scheme. To do so, 

• in time : the predictor- corrector scheme is used for the second order of accuracy in 
time, i.e. the numerical error is 0{dt^), instead of 0{dt^) 

• in space : in the same way, to achieve the second order scheme the linear distributions 
of quantities (conservative variables) inside cells are used. 

— Multi- dimension problem. At each timestep, we compute and sum the fluxes for each cell 
with all its neighbours and update the value of conservative variables with the total flux. 
Thus, we do not use operator splitting (dimensional splitting) and thus avoid the numerical 
artifacts introduced by this method, e.g. artificial spatial asymmetry. 

— Grid boundaries. To treat grid boundaries, we use the method of ghost cells. We 
include the two additional cells on either end of grid in each direction, and set the quantities 
in these cells at the beginning of each timestep. For simplicity, we set the quantities in 
ghost cells to be equal to these in the nearest "real" cell, thus implementing nonreflecting 
boundary conditions (outflow boundary). This physically correspond to boundary which 
does not reflect any wave, which is consistent with expansion into vacuum. 

— Vacuum treatment. In our simulations we deal with spatially finite systems expanding 
into vacuum. Thus computational grid in the Eulerian algorithm must initially contain both 
system and surrounding vacuum. To account for a finite velocity of expansion into vacuum, 
which equals c for infinitesimal slice of matter on the boundary, we introduce additional 
(floating-point) variables in each cell which keep the extent of matter expansion within a 
cell, having value 1 for the complete cell, for a cell with vacuum only. The matter is 
allowed to expand in the next vacuum cell only if the current cell is filled with the matter. 
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VI. RESULTS AND DISCUSSION 



In the article we apply the hydrokinetic model for an analysis of the space-time picture of 
Au+Au collisions at the top RHIC energies. Such an analysis provided in the evolutionary 
models of heavy ion collisions have to be based on a detailed description of the pion and kaon 
femtoscopic scales and must also describe well the absolute values of the spectra (not only 
spectra slopes) of the particles. As it was noted in Ref. jo], the following factors favor the 
simultaneous description of the mentioned data: a relatively hard EoS (crossover transition 
between hadronic and quark-gluon matters, not the first order phase transition), the pre- 
thermal transverse flows developed to thermalization time, an account for an "additional 



portion" of the transverse flows due to the shear viscosity effect [10|, a correct description 
of the gradual decay of the system at the late stage of the expansion. All these factors are 
included in the presented version of the HKM. 

We use both the Glauber-like (Section IIB) and CGC-like (Section IIC) initial conditions. 
In the former case the mean transverse radii, defined by ([1]) is Rt = 4.137 fm for the top 
RHIC energy. The best fit for the Glauber IC is reached at the following values of the two 
fitting parameters related to the proper time r = 1 fm/c: eo = 16.5 GeV/fm^ ((e) = 11.69 
GeV/fm^) and parameter of the initial velocity defined by ([1]), a = 0.248 {{vt) = 0.224). 
In the case of the CGC-like initial conditions = 3.88 fm, the fitting parameters leading 
to the best data description are eo = 19.5 GeV/fm^ ((e) = 13.22 GeV/fm^) and a = 0.23 
{{vt) = 0.208). The parameters a for the initial transverse fiows are somewhat larger than 



17| . The reason 



they are for the free streaming approximation of the pre-thermal stage 
is, as it is explained in Section II, that the fitting parameter a is related to the "unknown 
portions" of fiows, caused by the two factors: a developing of the pre-thermal fiows and the 
viscosity effects in the quark-gluon plasma. In addition, an account of the event-by-event 
fiuctuations of the initial conditions also leads to an increase of the "effective" transverse 
flows, obtained by averaging at the flnal sta ge, a s compare to the results based on the initial 



conditions averaged over initial fluctuations 



4l| . Since we use the later kind of IC, it should 



lead also to an increase of the effective parameter a. 

As it was discussed in Sections III, the chemically non-equilibrated evolution at the late 
stage, T < Tch = 165 MeV, is not characterized by a simple EoS, like p = p{e,fiB), in 
our calculations the pressure in this domain depends on 360 variables: energy density and 
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e [GeV/fm^] 



FIG. 1: (Color online) Equation of state p{e) used in the HKM calculations. The solid black line 
is related to the chemically equilibrated phase, taken from lattice QCD results as described in 
Sec. IIIA, while grey region consists of set of the points corresponding to the different hadron gas 
compositions at each e occurring during the late non-equilibrium stage of the evolution. The dashed 
line denotes EoS for the chemically equilibrated hadron gas and dotted line for the chemically frozen 
one, they are shown for a comparison. 

particle concentrations. In Fig. 1 we demonstrate the "effective" EoS at the temperatures 
around and below Tch- The points related to the later region characterize all diapason which 
the pressure gains at each energy density when the system evolves with the Glauber IC 
fixed above. We see that the pressure is differ from the "limited" cases: the chemically 
equilibrated and completely chemically frozen evolution (when the numbers of all (quasi) 
stable particles and resonances are conserved). At relatively large energy densities in a 
dominated space-time region the non-equilibrium EoS is harder than even in the chemically 
equilibrated case. This could reduce the out- to side- ratio for transverse interferometry 
radii. 

The results of the HKM for the pion and kaon spectra, interferometry radii and Rout/Rside 
ratio are presented in Fig. 2. Since the temperature and baryonic chemical potential at 
chemical freeze-out, which are taken from the analysis of the particle number ratios ^|, are 
is more suitable for the STAR experiment, the HKM results for kaon spectra are good for the 
STAR data but not so much for the PHENIX ones. Note also that, in spite of other studies 

n 

(e.g., |4|), we compare our results for the interferometry radii within the whole measured 
interval of pt covered at the top RHIC energy. Finally, one can conclude from Fig. 2 that 
the description of pion and kaon spectra and space-time scales is quite good for both IC, 
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the Glauber and CGC It is worth noting, however, that the two fitting parameters a and 
eo are various by 10-20% for different IC, as it is described above. 

The special attention acquires a good description of the pion and kaon longitudinal radii 
altogether with Rout/ R side ratio, practically, within the experimental errors. Such an achieve- 
ment means that the HKM catches the main features of the matter evolution in A+A colli- 
sions and correctly reproduces the homogeneity lengths in the different parts of the system 
which are directly related to the interferometry radii at the different momenta of the pairs 
{2!. In this connection it is valuable to show the structure of the emission function for pions 
and kaons. 

In Fig. 3 we demonstrate the space-time structure of the particle emission at the Glauber 
IC for different transverse momenta of particles, longitudinal momenta is close to zero. The 
space-time picture of particle liberation is quite different for different transverse momenta: 
for the soft particles the maximal emission occurs close to the cental part and happens at 
relatively later times, while the most of the hard particles are emitted from the periphery 
of the system at early times. In fact (see also IgI. I42I). the temperatures in the regions of 
the maximal emission are quite different for different pti they are for pions: T ?^ 75 — 110 
MeV for pr = 0.2 GeV/c and T ^ 130 - 135 MeV for pr = 1.2 GeV/c. So, if one uses 
the generalized Cooper-Frye prescription j^, 42| applied to the hypersurfaces of the maximal 
emission, these hypersurfaces will be different for the different particle momenta and does 
not correspond to common isotherm jo, I42I . 

One can see in Fig. 3, the top plots, that at equal transverse momentum pt the maximal 
emission of kaons happens earlier than pions as one can expect since the kaons interact 



weaker. At the same time the kaon interferometry radii in 



'ig. 2 follow approximately to 

the pion radii, demonstrating the approximate mT-scaling |43| with deviations to the slightly 
bigger values than pion radii have. The explanations can be gain from the middle row in Fig. 
3 where the comparison is done for the same transverse mass of pions and kaons. Then the 
maxima of pion and kaon emissions become closer and the majority of kaons leave system 
even somewhat later than pions at the same m-r, opposite to the comparison at the same 
Pt- Since in simplest situations the homogeneity lengths for bosons depend on ttit 
could say that the approximate mj'-scaling could indicate the similarity of the freeze-out 
picture for kaons and pions. However, probably, such a conclusion is very approximate since 
the real structure of the emission processes in A+A collisions is quite complicated as one 



43] . one 
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FIG. 2: (Color online) The transverse momentum spectra of (a) negative pions and (b) negative 
kaons, all calculated in the HKM model. The comparison only with the STAR data are presented in 
the separate small plots. The interferometry radii: (c) Rout, (d) Rside, (f) R-long and (e) Rout/Rside 
ratio for 7r~7r~ pairs and mixture of K~K~ and K'^K^ pairs. The experimental data are taken 
from the STAR 
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45| and PHENIX 



46|-|48| Collaborations. 



can see from the details in Fig.3. 
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't [«m] [tin] 

FIG. 3: (Color online) The (/)p-integrated emission functions of (a,b,c) negative pions and (d,e,f) 
negative kaons with different momenta: (a,d) pj- = 0.2 GeV, (b) px = 0.85GeV , (e) px = 0.7GeV, 
(c,f) pt = 1.2GeV at the Glauber IC. The values px in the middle row (b,e) correspond to the 
same transverse mass for pions and kaons uit = 0.86 GeV. 

VII. CONCLUSIONS 

The hydro-kinetic model is developed for a detailed study of the matter evolu- 

tion and space-time picture of hadronic emission from rapidly expanding fireballs in A+A 
collisions. The model allows one to describe the evolution of the QGP as well as the grad- 
ually decoupling hadronic fluid - a chemically non-equilibrium matter, where the equation 
of state is defined at each space-time point and accounts for decays of resonances into the 
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non-equilibrated medium. 

The HKM is applied to restore the initial conditions and space-time picture of the matter 
evolution in central Au+Au collisions at the top RHIC energy. The analysis, which is based 
on a detailed reproduction of the pion and kaon momentum spectra and measured femto- 
scopic scales, demonstrates that basically the pictures of the matter evolution and particle 
emission are similar at both Glauber and CGC initial conditions (IC) with, however, the dif- 
ferent initial maximal energy densities: it is about 20% more for the CGC initial conditions. 
The initial pre-thermal flow is slightly less for the CGC IC. The main factors, which allows 
one to describe well simultaneously the spectra and femtoscopic scales are: a relatively hard 
EoS (crossover transition and chemically non-equilibrium composition of hadronic matter), 
pre-thermal transverse flows developed to thermalization time, an account for an "additional 
portion" of the transverse flows due to the shear viscosity effect and fluctuation of initial 
conditions, a correct description of a gradual decay of the non-equilibrium fluid at the late 
stage of expansion. Then one does not require the too early thermalization time, < 1 
fm/c, to describe the data well. All these factors are included in the presented version of 
the HKM and it allows one to describe observables with only the two parameters. 

An analysis of the emission function at the top RHIC energies demonstrates that the 
process of decoupling of the fireballs created in Au+Au collision lasts from about 8 to 20 
fm/c, more than the half of fireball's total lifetime. The temperatures in the regions of 
the maximal emission are different at the different transverse momenta of emitting pions: 
T ^ 75 - 110 MeV for pr = 0.2 GeV/c and T ^ 130 - 135 MeV for pt = 1.2 GeV/c. A 
comparison of the pion and kaon emissions at the same transverse mass demonstrates the 
similarity of the positions of emission maxima, that could point out to the reason for an 
approximate scaling. 

Summary: the advanced HKM tool allows one to describe the process of the fireball 
evolution and gradual particle liberation in agreement with underlying kinetic equations. 
Further developments of the hydrokinetic approach and an analysis of the data in non- 
central A+A collisions will be the subject of a follow-up work. 
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